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Abstract. We address joint photodetection as a method to discriminate between the 
classical correlations of a thermal beam divided by a beam splitter and the quantum 
entanglement of a twin-beam obtained by parametric downconversion. We show 
that for intense beams of light the detection of the difference photocurrent may be 
used, in principle, in order to reveal entanglement, while the simple measurement 
of the correlation coefficient is not sufficient. We have experimentally measured 
the correlation coefficient and the variance of the difference photocurrent on several 
classical and quantum states. Results arc in good agreement with theoretical 
predictions taking into account the extra noise in the generated fields that is due 
to the pump-laser fluctuations. 



1. Introduction 



Entanglement is a crucial resource in quantum information processing, quantum 
communication and quantum measurements. Indeed quantum correlations lead to 
important novel effects not achievable by using classically-correlated states, i.e. states 
characterized by correlations that may be established by using local operations and 
classical communication. Quantum information has been initially developed for discrete 
quantum variables, i.e. quantum bits, which can be implemented optically by means of 
polarization single-photon states. However, much attention has been recently devoted 
to continuous variable (CV) regime and to multiphoton states of light. Continuous- 
spectrum quantum variables may be easier to manipulate compared to quantum bits 
by means of linear optical circuits and homodyne detection P3 |2l E] : this is the case of 
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Gaussian states of light, e.g. squeezed- and twin-beam. By using CV one may carry 
out nonlocality experiments [3], quantum teleportation [S] and generation of multimode 
entanglement The concepts of quantum cloning 17] and entanglement purification 
[S] have also been extended to CV, and secure quantum communication protocols have 
been proposed jH]- 

Ideal features for implementing quantum information experiments are the 
availability of bright and stable entanglement sources, based on degenerate or 
nondegenerate optical parametric processes, and the possibility of an effective 
characterization of entanglement. In the case of CV Gaussian entanglement quantum 
correlations may be discriminated from classical correlations by using homodyne 
detection. However, homodyne detection requires an appropriate mode matching of 
the signals with a local oscillator at a beam splitter, a task that may be particularly 
challenging in the case of pulsed optical fields. On/off photodetection may be also used 
to characterize Gaussian states, but its use is limited to states with a small number of 
photons HUE]. 

For the reasons of above, in this paper we assess the use of intensity measurements, 
in particular joint photodetection, as a method to discriminate classical correlations from 
entanglement A simple intensity-based measurement, including the measurement 
of difference photocurrent, cannot provide a complete characterization of entanglement. 
However, we show that for intense beams of light the detection of the difference 
photocurrent may be used, in principle, in order to reveal entanglement, while the 
simple measurement of the correlation coefficient is not sufficient. In particular, joint 
photodetection can be useful to discriminate the entanglement of twin-beam from 
correlations of thermal sources in the mesoscopic regime. We have experimentally 
measured the correlation coefficient and the variance of the difference photocurrent 
on several classical and quantum states. Results are in good agreement with theoretical 
predictions if one takes into account the extra noise in the generated states caused by 
the pump laser fluctuations. 

This work may also contribute to the recent debate on the use of classical and 
quantum correlations in imaging and on the necessity of entanglement for extracting 
the information ^Hl- Our results indicate that any method only based on correlation 
measurements cannot be improved using entanglement instead of classical correlations. 

The paper is structured as follows. In Section El we theoretically analyze the joint 
photodetection of classically and quantum correlated fields. In section 01 we present the 
experimental results obtained for a quantum (twin beam) and a classical (thermal light 
divided by a beam splitter) light. In Section HI we discuss the experimental results and 
draw our conclusions in Section El 

2. Quantum versus classical correlations 

Our aim is to assess the use of joint photodetection as a method to discriminate classical 
correlations from entanglement. The scheme we are going to consider is the following: 
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two modes of radiation, say hi and a 2 , are independently measured by two photodiodes, 
and the resulting photocurrents ih\ and m 2 are then electronically manipulated and 
analyzed. In the following we first investigate the use of the correlation function as an 
entanglement marker, and then pass to consider the difference photocurrent, of which 
we analyze both the variance and the distribution as a whole. The different markers are 
compared in order to discriminate entangled twin-beam (TWB) of radiation from i) a 
two-mode factorized coherent state showing no correlations, and ii) a two-mode thermal 
beam, showing classical correlations only. Since entanglement of TWB is a monotone 
function of its energy the comparison are performed for fixed mean number of photons 
of the involved signals. 

Entangled twin-beam Rx = \X))((X\ are obtained in quantum optics from 
(spontaneous) parametric downconversion (SPDC) in second order nonlinear crystals. 
The expression in the number basis is given by 

\x)) = vTh^fxx i*) 1 ® \ k h , (i) 

k 

where \k)j denotes a Fock number state in the Hilbert space of the j-th mode. The 
parameter x satisfies \x\ < 1 and may be taken as real without loss of generality. The 
value of x depends on the crystal length and on the nonlinear susceptibility, whereas the 
mean photon number of the TWB is given by ({X\rii + n 2 \X)) = 2N, where hj = ajaj 
with j = 1,2, and N = x 2 /(l — x 2 ) is the mean photon number of each beam. As 
a benchmark for uncorrelated classical signals we consider a two-mode coherent state 
of the same energy of the TWB, i.e. R a = \a)u(a\ <S> |«)22(«| with \a\ 2 = N. On 
the other hand, as a reference for classically correlated signals we consider the state 
obtained by sending a thermal state on a balanced beam-splitter whose second port is 
left unexcited. In general, if we mix a quantum state g with the vacuum in a beam- 
splitter of transmissivity r, the outgoing state is described by the density matrix 



R = ^2 T '' 1 i 1 -^ 2 " Qp+s,t+q 



stpq 




\s)(t\®\p)(q\ (2) 



where g^k = (h\g\h) are the matrix elements of the input state. In our case r = 1/2 
and the input state is a thermal state with 2N mean photon number, i.e g = v with 
Vh,k = <W(! + 2iV)~ 1 [2iV/(l + 2N)} k . We will denote the state obtained in this way as 
R v . As it can be easily seen by evaluating the eigenvalues of the partial transpose R 6 U , 
the state exiting a beam-splitter fed by a thermal state is never entangled, though it 
may show a high degree of classical correlations. 

We assume that photodetection is performed with quantum efficiency 77 and no dark 
counts. The probability operator-valued measure (POVM) of each detector, describing 
the statistics of detected photons, is thus given by a Bernoullian convolution of the ideal 
number operator spectral measure P nj = \nj){rij\ 



m 3 
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r n mi <g> n m2 



with j = 1,2. The joint distribution of detected photons p(m l5 m 2 ) can be evaluated by 
tracing over the density matrix of the two modes, i.e. p(mi,rri2) = Tr 

while the moments {m\rn\) = Tr Rrn\m\ 
of the operators 



v 

m 3 



of the distribution are evaluated by means 

£ mj n mj = (i - G ^ K) . ( 4 ) 



n,=0 



where 




(5) 



Of course, since they are operatorial moments of a POVM, we have, in general, mP- ^ rhP. 
The first two moments correspond to the operators 

rhj = rijhj 



rn 



2 _ ^2-2 



77-n- +7^(1 -Vj)nj ■ (6) 



3 '3 3 



As a consequence, the variances of the two photocurrents are larger than the 
corresponding photon number variances. We have 

a 2 ( mj ) = (m 2 ) - (rh,) 2 = a 2 ( nj ) + Vj {l - -q^hj) . (7) 

The correlation coefficient is defined as 

((mi - (mi))(m 2 - (m 2 ))) , . 

£ ~ 1 — w — \ y°) 

a{mi)a{m 2 j 

where rhj and a 2 {mj) are given in Eqs. (jUJ) and (J7J) respectively. Of course, for factorized 
coherent states we have e a = 0, while for the TWB and the thermal states we have 

£X v/(l + 7 7l JV)(l + 7 72 iV r ) ^ v/a + ry^a + rfeiV) 
which, for r/i = 772 reduce to 

_ (1 + AQ?? _ jvg nm 

£x " TTW £v ~ T+^n ■ (10) 

As it is apparent from Eqs. © and (jlUJI the correlation coefficient cannot provide 
a reliable discrimination of classical and quantum correlations for a mean number 
of photons larger than few units. As a consequence, any imaging system based on 
coincidence detection, cannot be improved by using entanglement. 

Let us now consider the quantity obtained by subtracting the two photocurrents 
from each other, i.e. the so-called difference photocurrent D = rh\ — rh^- The statistics 
of the outcome can be obtained as p(d) = Tr R 6^ where the POVM is given by 

ti q+d ®ti q d>0 

e d = }]{ fi g <s>n q d = o , (11) 
u q ®u q+d d<o 
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with n„ given in Eq. (J3J). The moments of the distribution can be obtained from the 
operators 

D = ^2dQ d = rjifii - rj 2 h 2 , 



D 2 = ^ d 2 <d d = fani - rj 2 h 2 ) 2 + ^(1 - rji)hi + r} 2 {l - n 2 )n 2 



(12) 
(13) 



which also provide the variance of the difference photocurrent cr 2 (d) = (D 2 ) — (D) 2 . 
For the class of states under investigation the difference photocurrent is distributed as 
follows 



p Q (d) = e -(*+*)* I w (2Ny/TMd 

n=0q=n+\d\ V 7 \ 



1 q=n+\d\ 




Xd 



(14) 



(15) 



p v {d) 



1 



1 + 2N 



E 



111*12 



E 



q.q' 



N 



1 + 2N 



g+q' 




n=0 

x (i- Vl y(i- m y' j^, (i6) 

where I n (x) denotes a modified Bessel' function of the first kind, and the J quantities 
are given by 



J, 



ad 



J, 



vd 




(17) 



(18) 



d < 



n. 



or d > and over 



In Equation (JTfij) the sum are over q = n + |c?|,..., q' 
q' = n + \d\, q = n, .. otherwise. The distributions are symmetric for 771 = t] 2 and 
asymmetric otherwise. In Fig. ^we display the distributions p a {d), px{d), and p v {d) 
for different values of the parameters rji, r] 2 and N. As it is apparent from the plots, 
the distributions for a thermal or a coherent state are broader than for the TWB, as far 
as the quantum efficiencies are close to each other and their value is not too small. In 
order to quantify this statement more explicitly we have evaluated, by using Eqs. (JEJ) 
and (fl3*jl. the variance of the difference photocurrent for the three types of states. We 
have 

=(vi + V2)N r '^ 2 2 V N (19) 

al{d) = ( Vl - V2 ) 2 N 2 + ( Vl + V 2)N ^ 2 2 V N (20) 

... (21) 



a 2 x (d) = (Vi - V2YN 2 + (77! + n 2 - 2 VlV2 )N ^ 2 V (1 - rj)N 
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For rji = 7] 2 = 7] the variances for the two classical states are equal, and larger than for 
the TWB state: the difference being more pronounced the greater is the rj value. On 
the other hand, if the two quantum efficiencies are different, we have o~ 2 (d) < a 2 (d) and 
o\ < a 2 (d) for any value of the mean photon number N, whereas o~ 2 x (d) < cr^,(d) only 
for numbers of photon below the threshold value 

»* = • < 22 > 

vrn - V2) 2 

In other words, for equal quantum efficiencies the variance of the difference photocurrent 
is a good marker to discriminate between quantum and classical correlations, whereas 
for different quantum efficiencies this statement is true only for signals with a small 
number of photons. In Figs. EJwe report the variances a 2 (d) as a function of the mean 
number of photons for both r]i = 772 and 771 7^ 77 2 , whereas in Fig. El we show a 2 (d)/N 
for rji = r]2 = 1] as a function of 77. 

Let us now consider a situation in which the two beams under investigation contain 
more than two, say 2/^, modes of the field, while the correlations to be discriminated are 
still pairwise. This is a common situation in pulsed experiments where several temporal 
modes are simultaneously matched in SPDC, and are present in thermal beams as well. 
We assume that the modes are equally populated. The statistics of counts for each 

detector is described by a multimode POVM of the form 

00 

m s =0 s 

where lt m is the single- mode POVM reported in Eq. (JHJ). The statistics of the difference 
photocurrent between the two detectors is described by a 2/z-mode POVM of the form 
with fl n replaced by Q n . 

Since the modes entering each detector are independent on each other we have 
(rhj) — ► (Y.s^js) = (J-irrij) and o 2 (m j ) -> ]^ s er 2 (m js ) = ficr 2 ^), j = 1,2. As a 
consequence the expressions of the correlation coefficients are still given by Eqs. © 
with iV that should be meant as the total mean number of photons of the \x modes. 
As concerns the distribution of the difference photocurrent we have, in terms of the 
probability density 



s q a ,r s 



5(J2q s -n-d)5(J2r s -n)9(d) 



+ 5(^2qs ~ n)5(S2r s - n - d)0(-d) 



(24) 



where 6(x) is the Heaviside step function. Notice that in writing Eq. 
we have already used the fact that the correlations are pairwise i.e. that 
p(lii r i) Q2, r 2, Qfi, r n) = n s p(g s ,r s ). By exploiting the delta functions in we 
may write 

00 d+n d+n-ni d+n—qi—...—q fl -i 

i*o = E££" E 

n=0 <ji=0 <J2=0 9/j =0 
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n n—ri n— n— r^_i 

Y Yl P(<l^ r i)p( ( l2,r2)...p(q ll ,r fl ) (25) 

r 1= T2=0 »>=0 

for d > and an analogue expression (with q s <-> r s ) for d < 0. 
3. Experimentals 

We verified the validity of the theoretical analysis on both quantum and classically 
correlated light. 

3.1. Twin Beam 

The quantum state of light we consider is a pulsed twin-beam generated by a traveling- 
wave amplifier in non-degenerate configuration. The layout of the experiment is depicted 
in Fig. |U As the pump source we use a frequency-tripled continuous-wave mode- 
locked Nd:YLF laser regeneratively amplified at a repetition rate of 500 Hz (High 
Q Laser Production, Hohenems, Austria). The laser delivers ~7.7 ps pulses at the 
fundamental frequency and ~4.5 ps pulses at the third harmonics. We obtain intense 
spontaneous parametric generation in broadly tunable cones by injecting the pump field 
(X p = 349 nm) into an uncoated /3-BaB 2 4 crystal (BBO, Fujian Castech Crystals, 
Fuzhou, China) cut for type I interaction (cut angle: 34 deg) having 10 x 10 mm 2 cross- 
section and 4 mm thickness. The pump beam, which emerges from the laser slightly 
divergent, is focused by lens f\ of 50 cm focal length. The crystal tuning angle is 
33.1 deg and the visible portion of the cones projected on a screen beyond the BBO is 
shown in the inset of Fig. EJ We operate in a dichromatic configuration by choosing the 
frequency of the laser second harmonics (Ai = 523 nm) for the signal and consequently 
the frequency of the laser fundamental (A 2 = 1047 nm) for the idler. For alignment 
purposes, a portion of the fundamental beam emerging from the laser is injected in 
the crystal together with the pump beam so as to obtain a well recognizable spot of 
amplified seeded down conversion. The selection of the two components of the twin 
beam is performed by means of two pin-holes, P\ and P 2 , having suitable dimensions, 
located on the outputs of the seeded process. In order to decide the dimensions of the 
pin-holes, such to collect a single coherence area at a time, we have to determine the 
dimensions of the coherence areas of the generated fields. In Fig. |5] (left) we show the 
single-shot picture of a portion of the signal cone taken with a digital camera (model 
Coolpix 990, Nikon, resolution 1024 x 768), in which we can clearly distinguish the 
presence of the coherence areas. In the right part (top) of the same figure we show a 
magnified single coherence area around Ai (green light) and (bottom) the intensity map 
of a typical coherence area taken with a CCD camera (model TM-6CN, Pulnix, operated 
at high- resolution). It is easy to demonstrate that the dimensions of the coherence areas 
in the idler beam (IR) corresponding to the measured signal beam scale according to the 
ratio of the involved wavelengths so that the dimensions for the idler are doubled with 
respect to the signal [H] . Accordingly, as shown in Fig. HI to select a single coherence 



Quantum and classical correlations of intense beams of light 



8 



area on signal and idler, we locate two pin-holes (diameter ~ 3.5 mm, on the signal and 
diameter ~ 7 mm, on the idler) at a distance of 72.5 cm from BBO. The light selected 
by the pin-holes is then focused with two lenses (/3 and f±, focal length 25 mm) on 
two p-i-n photodiodes (Si 85973-02 Hamamatsu, 1 ns time-response, 500 /zm diameter 
sensitive area on the green and InGaAs G8376-05, Hamamatsu, 5 ns time-response, 500 
/j,m diameter sensitive area on the IR) having nominal quantum efficiencies rji = 0.92 
and i]2 = 0.78 respectively. The current outputs of the photodiodes are integrated over a 
synchronous gate of suitable time duration (40 ns) by a boxcar averager that is operated 
as gated integrator in external trigger modality. The boxcar output is digitized by a 
13-bit converter (SR250, Stanford Research Systems, with 50 mV full-scale) and the 
counts stored in a PC based multi channel analyzer (MCA). The measurements are 
performed by inserting a variable filter (VF in the figure) in front of the photodiode 
detecting the signal, and by carefully adjusting it to balance the quantum efficiencies 
of the two detection branches of the setup. The interpretation of the output data must 
take into account the presence of cut-off filters, inserted to eliminate residual pump 
and all stray light, the overall quantum efficiency of the detection apparatus results to 
be rji ~ r]2 = 0.67. We verify the linearity of the boxcar integrators and measure the 
conversion coefficients («i = 6.7182 x 10~ 8 V and a 2 = 8.3043 x 10~ 8 V) by linking 
the voltage output of the digitizer to the number of electrons forming the photocurrent 
output pulse of the detectors at each laser shot. The relations among the statistics of 
the number of photons incident on the detector, p P h{n), the statistics of the number 
of detected photons, p e i(m), and the statistics of the output voltages of the acquisition 
apparatus, p ou t(v), are given by 



being a the measured conversion coefficient mentioned above and C a normalization 
coefficient. If we limit our analysis to the first two moments of the distributions, the 
experimental outputs are linked to Eq.s (0) and (J7J) by 



where for the sake of clarity we have defined a 2 el im) = o~ 2 (m) and cr 2 h (n) = a 2 in) [see 
Eq. (0)]. Note that in general the statistical distribution for the measured outputs is 
different from that of the incident photons. However, in both our cases (quantum and 
classical), the statistical distributions of the detected photons and of the voltage outputs 
are thermal ones. 

In Fig. IH1 we show the recorded signal (left) and idler (right) outputs of the 
photodiodes as a function of the laser shot, together with the noise of the detectors. 
In Fig. we the corresponding normalized probability distributions are reported for 
the same data. By looking at the probability distributions in Fig. [7| we note that the 




(26) 



Pout(v) = Cp e i(am) , 



(27) 



V = aM = ar]N 



2 [v 2 o- 2 ph (n)+v(l-v)N] 



(28) 
(29) 
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statistics of the outputs are well fitted by multithermal distributions [T^j, that is the 
distributions obtained by the convolution of \x equally populated thermal modes 

_ exp (-vfi/Vr) v^ 1 
,M " (//-I)! X JV^W ' { ) 

where Vt = aM? is the mean output corresponding to the overall detected photons 
mean value Mt- Equation (|HU|) holds in the high-intensity regime, which is the present 
experimental condition. In fact, by using the measured conversion coefficients on the 
detection arms of signal and idler we get Mi = 7.225 x 10 6 and M 2 = 7.212 x 10 6 
as the mean number of detected photons. As it is well known from the theory of 
photodetection [TH], the number of detected modes can be interpreted as the ratio of 
the time characteristic of the measurement (in our case the time duration of the pulse) 
and the coherence time characteristic of the field to be measured (in our case the inverse 
of the temporal bandwidth of the spontaneous parametric down conversion) ^H]- The 
continuous lines superimposed to the histograms of the experimental data in Fig.[7|show 
the convolution integrals, optimized for the number of temporal modes, of the theoretical 
distribution in Eq. ()30|) with the system impulse response evaluated from a measure in 
the absence of incident light. As expected, the signal and idler distributions are well 
fitted by multithermal distributions having the same number of modes (fi = 14). Note 
that the probability distributions for signal and idler are very similar to each other. In 
order to stress the correspondence between signal and idler, we plot the output of the 
idler as a function of that of the signal (see Inset in Fig.[HJ). To compare the experimental 
results with the theoretical predictions, we first of evaluate the correlation function of 
the photocurrents as 

m = Ef = i MAO-foi))Mfc + j)- (V2))/K ^ 

a(vi)a(v 2 ) 

where the average operations are taken over K (typically K = 30000) subsequent laser 
shots. For j = 0, Eq. (|3*T|) gives the correlation coefficient e 

£ = ((VI - (Vl)) (V2 - fo 2 ») ; (32) 

a(vi)a(v 2 ) 

which should be compared with the theoretical predictions of Eqs.flSj) and (|10|) . In Fig. [HI 
we show the correlation coefficient for the data of Figs. |H1 and the contributions of 
the noise of the apparatus, {i.e. the variance of the impulse response in Fig. Ej), are 
subtracted from the measured variances of the experimental data. We get e = 0.97, to 
be compared with a theoretical value of about 1. Note that subsequent shots results to 
be uncorrelated. 

As it has been shown in Sectional the distribution of the difference photocurrent is a 
relevant marker of entanglement. In Fig.Elwe plot the distribution of the difference of the 
photoelectrons detected on signal and idler, i.e. p(d) = p(m s — = p(vi/cti — v s /a s ). 
The distribution appears almost symmetrical and centered at zero, which indicates 
both accurate balance of the detectors' quantum efficiencies and high correlation in 
signal/idler photon numbers. The variance, as evaluated from the data, once the 
variance of the noise is subtracted, turns out to be o~\{d) = 2.124 x 10 11 . 
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3.2. Thermal Light 

To investigate joint photodetection for classically correlated light, we modify the 
experimental setup according to Fig. EH Pseudo-thermal light has been generated by 
inserting a moving ground-glass diffusing plate in the path of the second-harmonics 
output of the laser (A = 523 nm). A portion of diffused light is selected with an iris 
(in Fig. irU|) and then sent to a 50% cube beam splitter. The temporal statistics of the 
generated light can be described by the same statistics as in Eq. EDI [17J, in which the 
number of modes can be varied by changing the dimension of the iris in order to collect 
more than one spatial coherence area. The beams emerging from the beam splitter are 
then detected by the same apparatus used for the twin beam, where the pin photodiodes 
are now identical (model S3883-02, Hamamatsu, rj ~ 0.71, nominal) since the two beams 
are at the same frequency. The mean number of detected photons on the two beams 
are M 1 ~ M 2 ~ 2.22 x 10 8 . 

In Fig. ^2 we show the normalized probability distributions for the detected 
photons. The continuous lines superimposed to experimental data in Fig. [TT] are the 
best fits of the data obtained for 15 modes. As in the case of the twin beam, the 
two histograms are very similar and suggest a high degree of correlation that is easily 
verified evaluating the value of the correlation function. In the Inset of Fig. we plot 
the two voltage outputs of the beam splitter one versus the other, and in right part the 
correlation function for the classical beams in which again the contributions of the noise 
of the apparatus have been subtracted from the measured variances of the experimental 
data. We get e = 0.995 to be compared with a theoretical value of about 1. 

In Fig. H31 we plot the distribution of the difference of the photoelectrons detected 
on the two arms of the beam splitter. Again the distribution appears symmetrical and 
peaking at zero. The variance, as evaluated from the data upon subtraction of the noise, 
is al{d) = 4.097 x 10 13 . 

4. Discussion 

The experimental results discussed in SectionOJare obtained by keeping the values of the 
quantum efficiencies as close each other as possible. Therefore, they must be compared 
with the expected values for equal quantum efficiencies and with the shot-noise level 
for the intensities we are working at. The theoretical values are o~\{d) = 4.769 x 10 6 
and a 2 a (d) = 1.444 x 10 7 for the TWB and a 2 v {d) = a\ = 4.446 x 10 8 for the classically 
correlated thermal light. In order to obtain a realistic comparison between theory and 
experiment, we have to take into account the presence of noise that unavoidably affects 
the experimental data. We identify two main sources of noise: first of all, the difference 
between the overall quantum efficiencies on the two detection branches. In fact, although 
the experimental procedure was optimized so as to obtain the best balanced t] values, a 
small residual difference cannot be excluded, and, as we will see, a small balance error, 
even local across the beam to be measured, produces a relevant difference in the values 
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of a 2 (d). On the other hand, we have to take into account the unavoidable fluctuations 
of the laser source which affect all the fields under investigation. In fact, the pulsed 
pump field is not a plane wave having constant amplitude. Rather, its statistics is more 
realistically modeled by a Gaussian distribution, i.e. a Poissonian distribution affected 
by an excess noise [T£] 



p p (n) 



exp 



n 



(33) 



where a 2 



(n p ) + 5'noise an d & n i S e — % 2 i n p) 2 is the increase of the variance due to 



and 5 2 

fluctuations; the quantity x measures the amount of such a deviation. We will evaluate 
the influence on the generated beams of the excess noise in the pump by evaluating the 
error propagation. 



4-1. Imbalance of the quantum efficiencies 

To evaluate the modifications of the experimental results due to imbalance in the 
quantum efficiencies of the two branches, we equate the experimental results for 
a 2 (d) with the theoretical predictions for unbalanced quantum efficiencies of the 
photodetectors (see Eqs. (}2T)j) and (|2Tj) for 771 7^ r] 2 ). In the case of TWB, we obtain 
0.12 < 1 771 — 772 1 < 0.22 and in the case of classical field 0.05 < \r)i — 772 1 < 0.12. These 
values are too large to be reconciled with the high symmetry of the measured p(d) 
(see Fig. IT3J) . We can thus conclude that simply including a difference in the quantum 
efficiencies on the two detection branches is not sufficient to account for the experimental 
data. 



4-2. Fluctuations in the laser source 

We evaluate the influence of the excess noise of the third-harmonics pump pulse on the 
generated beams. 

Starting with the SPDC, we recall that the mean photon number in each component 
of the generated twin beam is given by 

N x = sinh 2 (ga p L) , (34) 

where g is a coupling constant, L is the interaction length inside the crystal and 
p = \/N p / (ApTp), being N p the mean photon number, A p the cross section and r p 



a 



the temporal duration of the pump pulse. By applying the error-propagation theory to 
Eq. ([34)1 . we get for the excess noise in the single mode of signal (idler): 

2 



5 x( n ) = °l 



fdN x y , 2 ^g 2 L 2 N 2 x 
— + x 2 J ^arcsinh 2 A/iVx 



x 2 ^arcsinh 2 v / iVx , (35) 
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where we used Eq. (|34j) and the final approximation holds for Af p >l. In the case of a 
mult it hernial beam composed by [i modes, Eq. ()35|) becomes: 



5 2 x {n) = — — x 2 arcsinh 2 y — — . (36) 
The variance of the difference photocurrent can thus be corrected as 



2 / j\ 2 

rfin U 77!// riSfj, y r}2ii 




o~ x (d) = cr X sp (d) 2^"X 2 arcsinh 2 4 / — 2^"X 2 arcsinh 2 ^/— - , (37) 



which is a function of the parameter x. We now evaluate the amount of laser fluctuations 
(i.e. the value of x) needed to reproduce the experimental data. To this aim, we equate 
Eq. (fHTj) to Eq. (J2H), modified to consider the presence of fi modes in the measured field 

M 2 M 
° 2 x(d) = (771 - r/ 2 ) 2 ^ + (m +772 - 2771%)-^ , (38) 

7/1,2^ Vl,2 

and study the dependence of x on the value of the overall quantum efficiencies on the two 
detected fields. Notice that, from the experimental point of view, we have two possible 
choices for the value of N appearing in the theoretical formula, namely N = Mj/rjj 
with j = 1,2 indicating either signal or idler, in our experimental conditions Mi ~ M 2 
and the two conditions give very similar results. Figure El displays the values of x 
as a function of 771 and 772 (left), and the corresponding values of the corrected 0x(<i) 
as calculated from Eq. (J37|) (right). The horizontal plane in on the right represents 
the shot-noise level of the measure as calculated from Eq. (|19|). Starting from data in 
Fig. we can draw two conclusions: On one hand, the experimental data corresponds 
to an amount of laser excess noise equal to x ~ 2.24%, which is compatible with the 
fluctuations of a pulsed laser. On the other hand, we have that at the intensities used in 
our experiments we cannot reliably discriminate the measured a x sp (d) from the shot- 
noise level. In fact, the right part of Fig. ITU shows that a slight indetermination in the 
quantum efficiencies may considerably increase the variance above the shot noise level. 
Note that the inclusion of an added noise does not imply a significant modification of 
the variance of the beams, as the total variance of signal/idler can be written as 

o\ = M (l + a; 2 arcsinh 2 ^^j . (39) 

As the correction to unity is less than 3%, the measured distributions are still well fitted 
by the expected multithermal distributions. 

As concerning the thermal light experiments, by applying the same strategy, we 
find that the excess noise can be written as 

6 2 u (n)=2x 2 N 2 , (40) 

which is again a function of the laser fluctuations x. Again we equate the value of the 
measured cr 2 (<i) corrected for the added noise 8 2 {n) and study the dependence of x on 
771 and 772. Figure ITB1 displays the values of x as a function of 77! and 772 (left), and the 
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corresponding values of the corrected o^{d) (right). The horizontal plane on the right 
represents the shot-noise level of the measure as calculated from Eq. (|19|). The highest 
values of the laser fluctuations, which are found for 771 ~ 772, is x ~ 3.6% at most. In 
contrast with the case of the TWB, from Fig. [T^lwe see that the values of al(d) are 
always above the horizontal plane representing the shot-noise level of the measure. 

In order to check the plausibility of the calculated values of x, we perform a stability 
measurements on the laser, by simultaneously detecting the second- and third-harmonics 
outputs of the laser with two photodiodes. In Fig. UHl we plot the measured values of 
x as a function of the third-harmonics energy in arbitrary units. The marked energy 
intervals in the plot indicate the operating range of the measurements discussed above. 
The obtained values of x are in agreement with those calculated. 

5. Conclusion 

Establishing the existence of entanglement and discriminating between classically 
and quantum correlated states in the high-intensity continuous-variable regime is a 
challenging task motivated by the need of characterizing the nature of the correlated 
light and of understanding the real resources needed to achieve the results in specific 
situations. We demonstrate that the characterization in terms of correlation functions 
is not satisfactory, as it gives similar results in both classical and quantum domain, 
whereas the measurement of the probability distribution for the difference photocurrent 
is in principle a good strategy. On the other hand, we demonstrate that in realistic 
high intensity conditions such a strategy cannot be reliably adopted, due to the 
unavoidable fluctuations of the laser source and slight imbalance of the detectors' 
quantum efficiencies. Indeed, by correcting the experimental data for these sources 
of noise, the data analysis leads to an agreement with the expected results. 

To achieve a more direct experimental demonstration we can follow two strategies. 
On one hand, we could work with identical quantum efficiencies, i.e. at frequency 
degeneracy, and use the same detection system on both parties of the correlated state. 
This could be done, for instance, by substituting the p-i-n photodiodes with a CCD 
camera. On the other hand, one may lower the intensity of the field to be measured, 
to decrease the sensitivity to the excess noise due to the pumping laser. Notice that, 
however, the possibility of lowering the intensity is limited by the amplifying capability 
of the electronic chain that manipulates the photodiode outputs. To overcome this 
limitation, one should switch to detectors with internal gain, such as photomultiplier 
tubes and hybrid photodetectors, taking into account that these detectors shows a 
low quantum efficiency of the photoelectric emission of the photocathodes which may 
compromises the overall visibility. 

In conclusion we have shown that difference photocurrent may be used, in principle, 
in order to reveal entanglement, while the simple measurement of the correlation 
coefficient is not sufficient. Our experimental results indicate that joint photodetection 
may be useful to discriminate the entanglement of twin-beam from correlations of 
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thermal sources in the mesoscopic regime. 
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Figure 1. Probability distributions p a (d) , px(d), andp„(d) for different values of the 
parameters iji, r\ 2 and N: the distributions for a thermal or a coherent state are broader 
than the corresponding distribution for the TWB, as far as the quantum efficiencies 
are close one each their and their value is not too small. 
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Figure 2. Variance <r 2 (d) of the difference photocurrent as a function of the mean 
photon number of the input signal. Left: for 771 = r)2 = 0.6; in this case o~ x {d) <C 
a a(d) = &i(d). Right: for ni — 0.5 and 772 = 0.7; for different 77's (T x {d) < (J 2 id) and 
a 2 a (d) < a 2 v {d) VJV, but a 2 x (d) < a 2 a (d) only for N < N th = 2 m m/(m ~ m? = 17.5. 




Figure 3. Ratio a 2 (d)/N between the variance of the difference photocurrent and the 
mean photon number of the signals as a function of the quantum efficiency, assumed 
to be equal for the two photodetectors. 




Figure 4. Experimental setup for measurements on the TWB: BBO, nonlinear crystal; 
/1-4, lenses; Fi_3, cut-off filters; VF, variable neutral filter; Pini_3, p-i-n photodiodes; 
P 1 2 j pin — holes; MCA-PC, multi-channel analyzer and data acquisition system. Inset: 
visible part of the downconversion cones. 
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Figure 5. Left: Single shot picture of a portion of the signal cone. Right-top: 
magnification of a single coherence area around Ai. Right-bottom: intensity map 
of a typical coherence area used to estimate its dimensions. 
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Figure 6. Left: Voltage outputs for the signal beam at Ai = 523 nm for a sequence 
of laser shots and noise. Right: Voltage outputs for the idler beam at A2 = 1047 nm 
for the same sequence of independently. 




Figure 7. Left: histogram of the intensity distribution of signal beam output at 
Ai = 523 together with the corresponding multithermal fit. Right: histogram of 
the intensity distribution of idler beam output at A2 = 1047 nm together with the 
corresponding multithermal fit . 
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Figure 8. Correlation coefficient between signal and idler beam as a function of the 
delay in the laser shot. Inset: values of the detected photons in the idler as a function 
of those in the signal in each laser shot. 




Figure 9. Experimental distribution of the difference photocurrent between signal 
and idler beams of a TWB. 




Figure 10. Experimental setup for measurements on classically correlated beams: 
BS, 50 % cube beam-splitter; D, moving diffuser; /1-4, lenses; VF, variable neutral 
filter; Pini_3, p-i-n photodiodes; P, iris; MCA-PC, multi-channel analyzer and data 
acquisition system. 
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Figure 11. Left:histogram of the intensity distribution of the transmitted beam 
together with the corresponding multithermal fit. Right histogram of the intensity 
distribution of the reflected beam together with the corresponding multithermal fit. 




Figure 12. Correlation coefficient between transmitted and reflected beams as a 
function of the delay in the laser shot. Inset: values of the detected photons in the 
reflected beam as a function of those in the transmitted beam at each laser shot. 




Figure 13. Experimental distribution of the difference photocurrent between the 
transmitted and the reflected beams of a pseudo-thermal beam impinging onto a beam 
splitter. 
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Figure 14. Laser fluctuations in experiments with TWB. Left: the amount of laser 
fluctuations a; as a function of the quantum efficiencies 771 and 772. Right: values of the 
corrected variance ~a\ (d) as a function of the quantum efficiencies rji and 772 ; the plane 
represents the shot-noise value. 





Figure 15. Laser fluctuations in experiments with thermal light. Left: the amount of 
laser fluctuations x as a function of the quantum efficiencies 771 and 772. Right: values 
of the corrected variance o^(d) as a function of the quantum efficiencies 771 and 772; the 
plane represents the shot-noise value. 
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Figure 16. Measured values of the laser fluctuations x for second- and third-harmonics 
outputs of the laser as a function of the third-harmonics energy in arbitrary units. 
The vertical lines delimitate the energy ranges of measurements performed on TWB 
an thermal light. 



